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Abstract 

The see-saw mechanism for neutrino masses based on the Grand Unification 
leads to the mass of the heaviest neutrino (~ u^) in the range (2 — 3) ■ 10"'^ eV and 
hence to a solution of the solar neutrino problem through the z/g — z/^ conversion. 
We suggest the existence of a light singlet fermion S which mixes predominantly 
with muon neutrino through the mixing mass m^^ ~ 0(1) eV. The introduction 
of such a singlet allows one (i) to solve the atmospheric neutrino problem via the 
T^ii ^ S oscillations, (ii) to explain the LSND result and (iii) to get two component 
hot dark matter in the Universe. The phenomenology of this scenario is considered 
and the origin of the singlet as well as of the scale m^g is discussed. 
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1 Introduction 



Grand Unified Theories containing heavy right handed neutrinos provide a natural framework 
for generation of small neutrino masses. Their generic properties are |I| 

(i) the (approximate) quark-lepton symmetry which relates the quarks mass matrices and 
the Dirac neutrino mass matrix at the Grand Unification (GU) scale lycu'- ~ ^g- This 
relation could be exact in case of the third generation so that mj^ = rrit at Aqu] 

(ii) the lepton number violation at the Grand Unification scale. 

This naturally realizes the see-saw mechanism of the neutrino mass generation with the 
Majorana masses of the right handed neutrinos at the GU scale: M/j ~ Aqu ~ 10^^ GeV. 
In this scenario, the mass of the heaviest active neutrino turns out to be in the range 
^ (1 — 3) ■ 10"'^ eV required for a solution of the solar neutrino problem via z/g — > Ur resonance 
conversion 0. Moreover, the z/g — 1/^ mixing can also be p, ^ in the correct range. The 
parameters of the z/g — mixing in this scenario can be such that the solar neutrinos undergo 
also long range vacuum oscillations on the way to the Earth ^. The interplay of conversion 
and oscillations can lead to observable effects [Q]. We will refer to this possibility as the Grand 
Unification scenario. 

The GU scenario does not allow one to explain the atmospheric neutrino problem. The 
existence of this problem and its oscillation solutions have been reinforced recently by high 
statistics observations of the zenith angle asymmetry and E/L dependence at the Super- 
Kamiokande as well as by the SOUDAN data 0. Furthermore, in the GU scenario the 
neutrinos are too light to explain the LSND [^] result and to be a hot dark matter relevant 
for structure formation in the Universe. 

We will show in this paper that a minimal extension of the GU scenario obtained by 
introducing only one light singlet fermion which mixes strongly with muon neutrino leads to 
solution of all these problems simultaneously. In what follows we will describe the scenario 
(sect. 2) and consider its phenomenology (sect. 3). In sect. 4 we will discuss the origin of 
singlet and its mixing. 

2 Scenario 

In analogy with the quark sector let us start with basis Ua = (z/g', z/^', z^r') where the leptonic 
mass matrices (both neutrino and charged lepton) have hierarchical structure with small mix- 
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ing. As in the GU scenario, we suggest that active neutrinos acquire a mass term via the 
see-saw mechanism, 

= -muM^^mJ . (1) 

Here mo and Mr are the Dirac and the Majorana mass matrices respectively. 

The scale in is identified with the grand unification scale A^^ and the Dirac mass of 
the third generation is identified with the top quark mass. This leads to rrirT ~ few x 10~^ 
eV. In case of the linear hierarchy of the mass eigenstates of Mji, the second diagonal mass 
m^^ ~ 10"^ eV and niee < m^/,- 

We suggest that on the top of the GU scenario there exists an additional light singlet state 
Us which couples predominantly to z/^/ in the basis Ua- 

mf,si^J^,Us + m^si^fus + h.c, (2) 

where 

m^s ~ 0(eV), russ < m^s- (3) 

The additional term m^^uj^iu^i can also be generated in some variants. We will take m^^ ^ 
rrins- The complete neutrino mass matrix can be parameterized as 



M,y = + m{i's) 



rriee 




rrier 
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rrier 




rrirT 















rriss 


J 



(4) 



In general, the elements rries and rrirs are non-zero. In this case the matrix can be reduced 
to the above form by additional rotation of usual neutrino components. We assume that this 
rotation is not substantially bigger than rotation corresponding to usual (CKM-) mixing. The 
mass rriss can be generated if z/^ couples to the right handed neutrinos: itIsrUs^r- This leads 
to rriss ~ ^IrI ^GU which gives m^s ~ lO^"*^ eV for nigR ~ 1 TeV. 

Let us introduce 4x4 matrix which diagonalizes (P, so that = U'^^Vi, where z/j 
{i = 1,2,3,4) are the neutrino mass eigenstates. Analogously, let f/' be a 4 x 4 matrix such 
that its 3 X 3 block diagonalizes the charged lepton mass matrix: la = U^a^a, = €.,fi,T) 
and f/^^ = Sas^as- Then the leptonic CKM matrix determining the fiavor mixing of neutrinos: 

= Uai^i equals 

Uai = U*[Uai, a = e,n,T,s i = 1,2,3,4; a = e', /x', r', s, (5) 
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The introduction of Ug with couphngs @ substantially changes the neutrino mass spectrum 
and lepton mixing of the GU scenario. 

(i) Since m^^ is much bigger than all see-saw contributions, u^r and Ug will form pseudo 
Dirac neutrino with mass ~ tti^^. The states u^r and Ug mix maximally in 1^2 and and the 
mass squared difference equals 

Aml4^ = ml- ml = 2m^s{m^^ + niss). (6) 

We fix Am|4 = (10~^ — 10~^) eV^, so that oscillations z/^ ^-^ Vg solve the atmospheric neutrino 
problem 0. For m^^ ~ eV, this leads to (m^/^ + mgg) ~ (0.005 — 0.05) eV. 

(ii) The introduction of Vg does not practically change the mass and mixing of the Ve' — ^r' 
sub-system. Indeed, diagonalization of the — Vg sub-system will give corrections to the 
mass terms of the z/g/ and Vr' suppressed by m^^/m^g ^ 10^'^ and m^^-jm^g <^ 1. Hence, 
^ai ~ ^ai = e',r','i = 1,3, where U^'^ diagonalizes the purely see-saw contribution, 
eq. ([^). Therefore as in the GU scenario [@ the v^, — conversion solves the solar neutrino 
problem 0. 

(iii) Unlike the i/g/ — i^^/ sub-system, the mixing between z/g/ — z/^/ as well as v^i — Vr' following 
from @ is suppressed in comparison with the see-saw prediction by factors : m^^/m^s ^ 10""^ 
and vTiT-T-jmy^g <^ 10~^ respectively. Therefore, the contribution of the neutrino mass matrix 
to lepton mixing in e — /i and /i — r sectors is negligibly small. The mixing will be determined 
mainly by the mass matrix of the charged leptons, 

can be chosen to explain the LSND result. 

Thus, we can write the matrix f/^^ which diagonalizes and therefore describes the 
contribution of the neutrino mass matrix to the lepton mixing as 
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(7) 



where sin0 = describes the z/g' — v^i mixing. For the charged lepton matrix ||f/^all 
assume a form typical for the quark mixing. Then the total lepton mixing matrix defined in 
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(I) can be written as 



U 



u;^,c-u;^,s u;^,/V2 u;^,s + u;,,c u;^,/V2 

u:,,c-u:,,s u:^,/V2 u:,,s + u:,,c u:^,/V2 

V -1/V2 J 




Here c = cos 6,s = sin 6. Notice that in our approximation has no admixture in the hght 
neutrino states z/i and z/3. Moreover, 



The resulting mass spectrum and mixing pattern are shown in Fig. 1. It is seen that this 
mass spectrum simultaneously solves the atmospheric neutrino problem, explains the LSND 
result and provides two component HDM in the Universe. 

3 Phenomenology 

For generic values of parameters the suggested 4 neutrino scheme has three different scales of 
the mass squared differences with the following hierarchy: 



Moreover, one mixing is almost maximal, whereas other mixings are small. This simplifies 
the task reducing dynamics in all particular cases to dynamics of two neutrino systems. Phe- 
nomenology of this scheme has been considers partly in ^ also. 

3.1 Solar neutrinos 

For solar neutrinos, the heavy components z/2 and z/4 "decouple" leading to the averaged 
oscillation effect, while the combination of the light states 



with cos^ (pe = Uei + ^e3 S^ts couverted to its orthogonal state. The survival probability can 
be readily written as 



Ua2 = UaA, = 6, /J,, T. 



(9) 



Am?n > Ami > Ami 




(10) 



P,y^ = COS^ 0eP2 + 77 Sin^ 0, 



(11) 
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where P2 is the two generation survival probabihty for z/^ (orthogonal state); here we have 
used U^2 = ^e4 = 0-5 sin^ 0e following from eq.(||). Strong bound on the angle 0e follows from 
the short range laboratory experiments. Indeed, for the short range experiments we find the 
survival probability 

P,^ = l- sin^ 20e sin^ (12) 



From BUGEY bound |Tg sin^ 20e < (4 - 8) ■ lO^^ for Amj^ ~ 0(eV) we get cos^ 0^ > 0.98 
and 0.5sin^(/)e < (0.5 — 2) ■ 10~^. Thus to a good approximation P^^ = P2{Aml^,sm'^ 29e), 
where sin^ 29e = 4:11^^11^^/ cos^ (pe- The corrections due to the presence of heavy states are 
very small. The electron neutrino survival probability does not significantly differ from the 
corresponding two active neutrino case which requires 

Am^3 ~ (4 - 10) X 10"^ eV^ ; AU^^U^^ ~ (3 - 10) x 10"^ (13) 
3.2 Atmospheric neutrinos 

The atmospheric neutrino deficit is explained mainly by the oscillations z/^ <-h> z/^. The transi- 
tion probabilities are small in other channels of oscillations. As follows from eq.(§), the muon 
neutrino has a small admixture of the lighter states ui and u^. This leads to oscillations due to 
large splitting Am^2 which get averaged out. Strong non averaged oscillations are stipulated 
by splitting Am24 while the splitting Am|]^ is irrelevant due to smallness of the z/g mixing in 
the heaviest states. As a result, we find the survival probability for the muon neutrinos: 

P., = l- cos' 0^ sin^ - 1 sin^ 20^ , (14) 



where cos^0^ = + ^^4) we used equality = f^^4- The last term in ([T^) is due 
to averaged oscillation contribution. The modifications of the usual two neutrino oscillation 
probability in (|T^) described by the parameter 0^ are strongly constrained by the short-range 
experiments. Indeed, for the the oscillations with Am1^ ~ Am2j ~ eV^ (i = 1, 3) we find the 
survival probability 

P,^ = l-sin2 20^sin2^^. (15) 



The negative results in the z/^ disappearance experiment by CDHS group imply |Tl| sin 20^ < 
0.1 for 1.0 eV^ < Aml^ < 10 eV^. Thus the last term in (|T^) is less than 0.05 and P^^ reduces 
approximately to two neutrino case. In particular, the double ratio in the atmospheric neutrino 
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fluxes is determined by the survival probability P^^. The analysis of the data gives relevant 
ranges of the neutrino parameters 

Ami, ~ (2 - 10) X 10-3 eV^ ; 4U%Ul, ~ (0.8 - 1) (16) 

which are similar to those for u^ — Ut-. 

The z/^ — Vr and — Vg oscillation solutions are rather similar as far as the low energy 
charged current data are concerned. There are crucial points which distinguish the solutions. 

1. The neutral current induced events are unaffected if the atmospheric neutrinos oscillate 
to Vr but are reduced in the case of the — Ug oscillations [0 . The neutral current reactions, 
in particular uN uNtc^, give main contribution to the sample of isolated events. Use- 
ful measure of this is provided by the ratios {i)N^o/Nf, {ii)NT,o /N^ and {ni)NT,o /Nmuiu-ring- 
Here N^^^^^^o denote the standard e-hke, /x-like events and number of events with isolated pions 
produced mainly through neutral current interaction. The Nmuiu-ring result mainly through 



the charged current induced reactions [13|. The ratios can differ by (30-50)% in two relevant 
cases if the survival probability is around 0.6 ||13|. It should be possible to see this differ- 
ence and distinguish these two possibilities with increased statistics in the Superkamiokande 
experiment. 

2. Events induced by the high energy neutrinos, like upward going muons (throughgoing 
and stopping) are affected by matter in the case of — Ug oscillations and there is no matter 
effect in z/^ — Ur channel. For z/^ — Ug oscillations the matter effect leads to peculiar zenith angle 
(6) distribution of the muons [0, ITSf. The distribution has two dips at cos© = —0.8 — 0.3 
and cos = —1.0 — 0.8. The latter dip is due to the parametric resonance in oscillations of 
neutrinos which cross the core of the Earth. In fact, MACRO, Baksan and SuperKamiokande 



data give some indications of such a distribution [|15|. 

3. Up-down asjTiimetry in neutral current induced events, i.e. uN uNtt^ can also 
provide a means to distinguish the two oscillation channels |]lj 



3.3 Short range — z/g and — oscillations 

The transition probability Pp^v^^ in the short range experiments which are sensitive only to the 
large mass splitting Am\^ Am^g ^ Am\^ ^ Am\^ is given by 

n • 2 n I -2 AlTi/i t 
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where according to (H) 

Positive signal reported by LSND implies 

sin^ 2<pe^, ~ (0.1 - 4) ■ 10"^ Am^^ ~ (0.3 - 2.5) eVl (17) 



Notice that 9^^ ~ 6*^ ~ Jnie/m^ gives for the effective mixing angle appearing in (|T 



sin^ 20e/, ~ 2 X 10"^ 

which is in the right range. 

For — Vj. short range oscillations we get similar two neutrino oscillation probability with 
Am^ = Am^4 and 

sin2 20^. = -4(f/,if/^i + Ur^U^^){Ur2U^2 + UrAU^A) ~ A{ul^,ul^, f . 

The depth of oscillations is determined basically by mixing of charged leptons in the selected 
basis. 

3.4 Supernova neutrinos 

Let us first find the level crossing scheme, i.e. the dependence of the eigenvalues of the 
Hamiltonian of the system on the density of the medium. The Hamiltonian can be written as 

M2 

H= + V, 

2E 

where M is the neutrino mass matrix in flavor basis and V = diag{Ve, V^, V^, 0) is the matrix 
of the potentials. 

The levels can be found in the following way. One diagonalizes first the strongly mixed 
heavy sub-system and z/^. This gives the levels i/g^? ^im- Then using smallness of all other 
mixings one gets the level crossings: z/g — 1^4^, — ^2my large densities and z/g — crossing 
at small density. There is no level crossing in the anti neutrino channels. 

In case of the complete adiabaticity one predicts the following transitions: 



<m ~ i^, + ^s)/V2, (18) 
J^f,^ ^T, (19) 

lyr^l^e- (20) 
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Thus at the surface of the star one will have hard (corresponding to the original Ur spectrum) 
electron neutrinos, u'^s with soft spectrum and the flux of sterile neutrinos which is about 1/12 
of the total neutrino flux. No changes in the antineutrino channels are expected. 

The production of the heavy elements due to r - processes in the supernovae imply that 



transitions u^- z/g, t'^ z/g are suppressed in the inner parts of the star [0]. As follows from 
the level crossing scheme the appearance of can be due to adiabatic transition — > z/g. 
The transition occurs above the region of the r-processes if Am^ ~ m^^ < 2 eV^. For larger 
mass square differences mixing angles of active neutrinos should be strongly suppressed. 

3.5 Primordial Nucleosynthesis 

The oscillations z/^ Vg and v^. lead to formation of the equilibrium concentration of the 
sterile neutrinos in the Early Universe. Therefore by the time of the primordial nucleosynthesis 
there were N^^ = 4 neutrino species. This is inconsistent with the bound Ni, < 2.6 based 



on low values of deuterium abundance |TU| but is allowed by recent conservative estimates 
N,y < 4.5 ||2y]. The appearance of sterile neutrinos can be suppressed if there is sufficiently 
large leptonic asymmetry in the Early Universe at the epoch with T < 10 — 20 MeV . However, 
this asymmetry cannot be produced in z/^ ^ z/^, z/^ — z/^ oscillations by mechanism suggested 
in because of the negative sign of the mass squared difference mg — in the scenario 



under consideration. 

4 Origin of scales 

Both the origin and lightness of sterile state can be understood in supersymmetric schemes 
p3| , p4| . The sterile neutrino may be identified with fermionic component of the Goldstone 
multiplet arising due to spontaneously broken global symmetry such as Peccei Quinn (PQ) 
symmetry or with some hidden sector field in string based models. Violation of the R 



symmetry in these theories could lead to coupling of the sterile state to ordinary neutrino 
The general discussion of the mass and couplings of such sterile state was given in |2^, |24 



Here we discuss a specific example which can reproduce the mass and mixing pattern discussed 
above. 

The example is based on the supersymmetric standard model with additional U{1)q sym- 
metry. We need to add singlet fields a, a' and y in order to implement the breaking of the 
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additional symmetry. The U{1)q charge is defined as 



Q = B- L-X . 

where B and L are the baryon and the lepton charges correspondingly. Keeping ultimate 
unification in mind, the X charge is assumed to commute with SO{10). The X charge can be 
normalized such that the common charge of the standard fermions is one unit. The X-charges 
of the Higgs fields {Hi, H2, a, a', y) are then chosen as (—2, —2, 2, —2, 0). Note that the U{1)q 
symmetry may be identified with the PQ symmetry. In addition to the standard Yukawa 
interactions, the following terms are allowed in the superpotential by the U{1)q. 

2 

W = \y{aa' - f) + M,N,N, + e.L.H^ + SH^H^^ . (21) 

IVIp 

where / and are the mass parameters, Mp is the Planck mass, and 6 = 0(1). The first term 
leads to the VEV (cr) = (cr') = / and thus the U{1)q breaking which generates the Goldstone 
multiplet S = (a — a' )/V2. The mass terms for the right handed neutrinos simultaneously 
break B — L and X but preserve U{1)q. It is natural to identify the masses Mj with the GU 
scale in the present context. The parameters allowed by the U{1)q automatically break the 
-R-parity. The last term generates the /i parameter of the required mag nitude: /i = 6p/Mp 
and also leads to coupling of S to neutrinos because of the presence of e^. We need to assume 
ei,3 £2 in order to ensure that S predominantly couple to u^. 

It is easily seen that the last two terms generate mixing of z/, S and Higgsinos. The block 



diagonalization of the mixing matrix gives |23] 



m^s = — ^ , (22) 



when heavy Higgsino states are decoupled. Here c ~ 0(1). 

The e2-term also generates mixing between i/^ and gaugino. This mixing in turn gives the 



contribution to the mass of the muon neutrino [^, m 



2 



rrif,^ = y^j mo (23) 

The above contribution dominates over the standard see-saw contribution to m^^ if the right 
handed neutrinos have masses at the GUT scale. The mass ttiq is completely determined by 
the parameters of the supersymmetric standard model if one starts with the universal scalar 
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masses at the GUT scale, mo is mainly determined in this case by the b Yukawa coupling hb 
and is typically given by 



msuSY \1Q7T^ M| 



Detailed analysis shows |26[ that mo ~ keV — MeV, when MSSM parameters are varied over 
its allowed range. 

The breaking of supersymmetry results in generation of the direct mass term rriss for 5*. 
The scale of rriss need not coincide with the susy breaking scale but could be much smaller, 



e. g. mgjjgy/^P Pll ''^susY / ^GUT depending upon the details of underlying theory. 



The niss is found ^ to be suppressed if {i) all the MSSM fields and some or all of the 
singlet fields have no scale kinetic terms or (m) if supersymmetry breaking is mediated by the 
gauge interactions and the MSSM fields as well as the singlets a, a', y and the right handed 
neutrinos do not directly couple to the sector in which supersymmetry is broken, [iii) Singlet 
S is identified with some of the modulie fields. The russ can arise radiatively in some cases 
p3| but its magnitude can be smaller than the neutrino-gaugino mixing contribution to m^^. 
The latter contribution then sets the scale for the atmospheric neutrino anomaly. 

The mixing mass m^s is in the eV range if €2/ f ~ 10^^^. The m^^ can be ~ 10~^eV 
provided 

{t^/^if ~ 10-^ - 10-s. (24) 

It follows then that 

£2 ~ (10^2 - 1) GeV ; / ~ (10^ - lO") GeV. (25) 

The value of / required here is included in the allowed range for the PQ symmetry breaking 
scale. The hierarchy 61^3 ^ t2 among the R breaking parameters can be explained by intro- 
ducing the additional (approximate) horizontal symmetry prescribing non-zero charges for the 
first and the third generation and zero charge for second generation and singlets. 

We have identified the singlet state with a supersymmetric partner of a Goldstone boson in 
the above example. In a more general context, singlet may remain light for other reasons. The 
desired value of the mixing mass can be obtained in such cases from the following relation: 

<X>pv m^/2Mpv 



A2 A2 



where < X >f 'i^s, the VEV of the field which breaks SUSY. This possibility can be realized if 
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e.g. the term 

^ (26) 

exists in the Kahler potential (it might be if unification with gravity occurs at the Grand 
Unification scale). The effective interaction ( p6|) can be generated due to exchange of the 
heavy superfield doublet which has the couplings in the superpotential: LS^u and XHd^u- 

5 Discussion and Conclusions 

The phenomenology of our scenario strongly differs from the phenomenology of the 4-neutrino 



schemes suggested in ^ |24|. There solar neutrinos are converted into sterile neutrinos, 
whereas the atmospheric neutrinos undergo — v^- oscillations. In contrast, in the present 
scheme the solar neutrinos are converted to active v.^ and atmospheric neutrino deficit is due 



to — Vs oscillations. 



The atmospheric neutrino problem gets resolved in terms of the z/^ oscillations 
also in the scenario proposed in |T^. However solar neutrinos are converted to the state 



which has a comparable admixture of the active and sterile components. These possibilities 
can be distinguished by studying the neutral current interactions in the solar neutrinos (SK, 
SNO experiments) as well as NC induced events (tt^) and zenith angle distributions in the 
atmospheric neutrinos. In fact in the model fl^ there is no explanation of the LSND result 
for the generic values of parameters. 

The hierarchy studied here was also proposed in purely phenomenologically. We have 
analysed this hierarchy in the context of specific see-saw models with Grand Unification. The 
difference between purely see-saw predictions and the present hierarchy can arise from the 
physics at the intermediate scale such as the PQ breaking. Also it is possible that all light 
scales are associated with GU scale and the scale of SUSY breaking. 

The same pattern of mixings and masses as proposed here was shown to arise in a model 
P7| with purely radiatively generated neutrino masses. The DM neutrino mass appears in 
one loop whereas the solar neutrino scale is generated in two loops. For certain choice of 
parameters the splitting between the two heavy states can solve the atmospheric neutrino 
problem. It is not always easy to include grand unification into such models which postulate 
an existence of new charged scalar bosons fc"^"*", 77+ and at the electroweak scale. This 
model can be distinguished by detection of the effects of these new bosons. 
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A see-saw type of scheme with the mass hierarchy similar to ours and based on the singular 
see-saw mechanism has been considered in also. The right handed neutrinos are required 



to be very light (~ keV) in [25| unlike the standard GU scenario pursued here. The eV mass 
scale for the neutrinos also does not come out naturally. 

In conclusion, the presence of a non-zero neutrino mass as hinted by the solar neutrino ob- 
servations may ultimately be pointing to the underlying grand unified structure. We proposed 
a minimal extension of this basic structure which also accommodates explanations for (i) the 
atmospheric neutrino anomalies (ii) LSND results and (iii) provides hot dark matter in the 
universe. The scenario suggested here can be tested in future solar neutrino experiments as 
well as in atmospheric neutrino observations at Superkamioka with increased statistics. The 
basic scenario proposed here could arise from an underlying supersymmetric theory. Thus 
verification of the scenario may indirectly probe supersymmetry as well as grand unification. 
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Fig. 1 Qualitative Pattern of neutrino masses and mixing explaining differ- 
ent neutrino anomalies. Boxes correspond to different mass eigenstates. The 
sizes of different regions inside each box determine the flavour content of the 
corresponding mass eigenstate. The electron (muon) flavours are represented 
by the weakly (strongly) hatched regions and black region represent the tau 
flavour. The sterile neutrino content is represented by white regions. Arrows 
connect eigenstates involved in solution of the solar (vq), atmospheric (ATM) 
and LSND problems. 
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